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A search for high redshift molecular absorption lines toward 
millimetre-loud, optically faint quasars 
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ABSTRACT 

We describe initial results of a search for redshifted molecular absorption toward four 
millimetre-loud, optically faint quasars. A wide frequency bandwidth of up to 23 GHz per 
quasar was scanned using the Swedish-ESO Sub-millimetre Telescope at La Silla. Using a 
search list of commonly detected molecules, we obtained nearly complete redshift coverage 
up to z a b s = 5. The sensitivity of our data is adequate to have revealed absorption systems 
with characteristics similar to those seen in the four known redshifted millimetre-band ab- 
sorption systems, but none were found. Our frequency-scan technique nevertheless demon- 
strates the value of wide-band correlator instruments for searches such as these. We suggest 
that a somewhat larger sample of similar observations should lead to the discovery of new 
millimetre-band absorption systems. 
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1 INTRODUCTION 

Millimetre-band (mm-band) molecular absorption systems along 
the line-of-sight to qu asars provide a powerful probe of cold gas 
in the early Universe. IWiklind & Combed ll994ii Il995l [l996b) 
have used molecular absorption lines to study a variety of prop- 
erties of the absorbers themselves (e.g. relative column densi- 
ties, kinetic and excitation temperatures, filling factors etc.). Be- 
sides information about the absorbers, important cosmological pa- 
rameters can be extracted from such data. Constraints on the 
cosmic microwave background temperature can be obtained by 
comparing the optical depths of different rotational transitions 
[e.g. CO(l-2), CO(2-3) etc.] (e.g. Wik lind & Combed Il996ch . 
Also, if the background quasar is gravitationally lensed, time de- 
lay studies can yield co nstraints on the Hubble constant (e.g. 
IWiklind & Combesl200ll) . However, these studies have so far been 
limited by the paucity of mm-band molecular absorbers. Only 4 
such syste ms are curre ntly known: z a bs — 0.685 toward TXS 
0218+357 
1413+135 
1504+377 
PKS 1830- 



Wiklind & Combes 1995) 



Wiklind & Combes 1997) 



0.247 toward PKS 
0.672 toward TXS 
Wiklind & Combedll996J) ari d z abs = 0.886 toward 
211 IWiklind & Combe Il998h . 
Quasar absorption lines can also be used to search for possi- 
ble variations in the fundamental constants. Detailed studies of the 
relative positions of heavy element optical transitions in 49 high 
redshift (0.5 < z a b s < 3.5) absorption systems favour a smaller 
fine structure constan t (a = e 2 /he) at t he 4. 1 a significance level 
jMurphv et al]l2001ct rWebb et alj|200ll) . The observed fractional 
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change in a [Aa/a = (-0.72 ± 0.18) x 10 -5 ] is very small 
and systematic errors have to be carefully considered. However, 
a thorough search for systematics has not revealed a simpler ex- 
planation of the optical results jMurphv et all200lj) . Independent 
constraints at similar redshifts are required and recent attention has 
focused on molecular absorption systems. 

Comparison of molecular rotational (i.e. mm-band) and cor- 
responding Hi 21 -cm absorption line frequencies has the poten- 
tial to constrain changes in a with a fractional precision ~ 10 -6 
per absorption system - an order of magnitude gain per absorption 
system over the purely optical methods. The ratio of the hyperfine 
(21 -cm) transition frequency to that of a molecular ro tational line 
is oc y = a 2 g p for g p the proton g-factor (Drinkwater et al. 1998). 
Thus, any variation in y will be observed as a difference in the ap- 
parent redsh i fts, A y/y « Az/(1 + z) = (z mo i - Zh)/(1 + Zmnl)- 
ICarilli et all feOOOl) and lMurphv et alj l2001bl) have obtained con- 
straints on Ay/y consistent with zero y-variation from spectra of 
PKS 1413+135 and TXS 0218+357. Currently, the major uncer- 
tainty in this mm/Hi comparison is that intrinsic velocity differ- 
ences between the mm and Hi absorption lines are introduced if 
the lines-of-sight to the mm and radio continuum emission regions 
of the quasar diff er, as is cer tainly the case for PKS 1413+135 and 
TXS 0218+357 iCarilli et alfcOOfj) . 

A statistical sample of mm/Hi comparisons is therefore re- 
quired to provide a reliable, independent check on the optical re- 
sults for Q-variation. One systematic approach to finding more mm- 
band molecular absorbers is to scan the frequency space toward 
a sample of millimetre-loud quasars. Indeed, the absorber toward 
PKS 18 30-211 was identified in this way by IWiklind & Combed 
ll996d) . With the assumption that molecular absorption will be as- 
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Table 1. The source list. Si. 4, S2.7,— £23 are the measured 1.4, 2 .7,... 230 GHz continuum flux densities in Jy. All values are obtained from 
the current version of th e Parkes catalogue jWright & 'otrupcek 19 9(1), except ^Phil lips et alj ll995l) . the 6 NVSS catalogue ICondon et alj|l99^) . 
iTerasranta et al J IT998 ), " Kova lev et alj JT*999l > and fromikoski. Lainela & Valtaoial l200Ct> . No optical magnitudes are published for these objects. 
Howe ver, note that the sight-lines to both B 0648—165 and B 0727— 115 suffer from moderate Galactic extinction, Ag jSchlegel, Fink beiner & Davis] 
1998). Also, B 0500+019 is known to have an intervening H I absorption system. Since we do not reliably detect continuum emission from this quasar 
(see Fig.0, we cannot place interesting limits on HCO+(l-2) or CO(2-3) absorption at 1 12.570 GHz and 218.227 GHz respectively. 



Quasar 




Coordinates (J2000) 


A B (mag) 








Radio flux densities (Jy) 










h m s 


d ' " 




SlA 


S2.7 


S5.0 


S'8.4 


Su 


S22 


Sse 


S230 


B 0500+019/J 05034 


203 


05 03 21.2 


02 03 05 


0.289 


2.10 


2.46 


2.04 


1.61 


1.36 d 


0.86 d 






B0648-165/J0650- 


1637 


06 50 24.6 


-16 37 40 


2.456 


1.70 


1.40 


1.02 


0.80 






0.9 a 




B 0727-1 15/J 0730- 


1141 


07 30 19.1 


-11 41 13 


1.271 


2.66 6 


1.95 


2.22 


3.36 


3.87 d 


2.91 d 


0.9 a 




B 1213-172/J 1215- 


1731 


12 15 46.8 


-17 31 45 


0.253 


1.50 


1.33 


1.28 


1.56 


2.97 d 


2.44 d 


1.5° 


ps 0.7 c < e 



sociated with significant optical extinction, one should select opti- 
cally faint quasars to increase the probability of detecting molecu- 
lar absorption. In this paper we present wide-band millimetre-wave 
spectra of the four millimetre-loud quasars which have not yet been 
optically identified: non-detections in the APM 1 and DSS 2 cata- 
logues imply tub > 20. In the following section we describe the 
observations and data reduction and present the wide-band spectra. 
In Section[3]we search for possible millimetre absorption systems 
close to our detection limits. We make our conclusions and discuss 
the future of our search technique in Section|4| 



2 QUASAR SPECTRA 
2.1 Observations 

We observed the quasars listed in Table Q in February 2002 with 
the 15-m SEST at La Silla, Chile. The receivers were tuned to 
single-sideband mode and typical system temperatures, on the T\- 
scale, were «250 K for the SESIS RX100 and RX150 receivers and 
ps340 K and 480 K for the IRAM RX1 15 and RX230 respectively. 
The backends were acousto-optic spectrometers with 1440 chan- 
nels and a channel width of 0.7 MHz. We used dual-beam switching 
with a throw of pa 12' in azimuth, and pointing errors were typically 
3" rms on each axis. 

We used the SESIS and IRAM RX230 receivers for the major- 
ity of our integrations, providing the advantage of a wide (1 GHz) 
bandwidth (the IRAM RX1 15 has a 0.5 GHz bandwidth). Since the 
backend response decreases sharply toward the edges of the band, 
we overlapped the bands by observing at intervals of 0.8 GHz to 
ensure uniform signal-to-noise ratio (S/N) in the final, combined 
spectra. The lowest frequency we could tune the RX100 receiver 
to was 78.8 GHz and, although the nominal range is 78-1 16 GHz, 
the backend configuration did not allow tuning to 81.2, 82.0 or 
82.8 GHz. Similarly, we could not tune the RX150 receiver to 141.7 
or 142.5 GHz. 



2.2 Data reduction 

For each 1 GHz- wide spectrum, we averaged typically 10 two- 
minute (reference-subtracted) scans of the source and subtracted 
a low order polynomial fit from the data to provide a flat contin- 
uum. To convert the data from the TJ-scale to the Jy-scale we fit- 
ted a third-order polynomial to the aperture efficiency of SEST as 
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Figure 1. The average flux densities, S v , for each 1 GHz- wide portion of 
spectrum. The error bars represent the rms derived from typically 10 scans 
at each frequency, v. Power-law fits with the given parameters were used 
to define the continuum flux to which the spectra in Figs. I2II5I have been 
normalized. 



a function of frequency: 

S v /Jy = 7.55 + 0.340(iVGHz) - 2.07 x 10" 
+5.32 x 10~ 6 (i7GHz) 3 . 



'(i//GHz 



(1) 



The intensity was calibrated using the chopper-wheel method 
which should be accurate to ±10% 3 . However, we found signifi- 
cant variations in the total flux density between each two-minute 
scan. Therefore, in Fig.Qwe present the mean flux density, S v , for 
each of the 1 GHz- wide spectra with error bars representing the rms 
from the contributing scans. 

For each channel, we generated a 1 a error from the rms in a 
window of width 27Y crr + 1 channels centred on that channel. We 
found that N CII = 20 provided reliable rms values but the results 
in Section l3~2l were insensitive to this parameter. In the 0.2 GHz 
regions where different 1 GHz-wide spectra overlapped, we re- 



1 Available at 

2 Available at 



http://www.ast.cam.ac.uk/~apmcat 



http://archive.stsci.edu/dss 



From the SEST handbook, http://www.ls.eso.org/lasilla/Telescopes/SEST 
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Frequency [GHz] 



Figure 2. Wide-band spectrum of B 0500+019. These data (and those in Figs. l3ll5l are normalized by the continuum flux of the quasar (taken as a power-law 
fit to the flux density data in Fig.0 and, for clarity, have been boxcar-smoothed over 3 channels. The mean S/N = 1.72 per (unsmoothed) channel. 
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Figure 3. Wide-band spectrum of B 0648—165 (see caption to Fig.|2|for details). The mean S/N = 5.73 per (unsmoothed) channel. 
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Figure 4. Wide-band spectrum of B 0727— 115 (see caption to Fig.|2|for details). The mean S/N = 9.63 per (unsmoothed) channel. 
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Frequency [GHz] 

Figure 5. Wide-band spectrum of B 1213—172 (see caption to Fig.|2|for details). The mean S/N = 3.40 per (unsmoothed) channel. 
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sampled the variance-weighted average using the largest channel 
width of the contributing spectra. Finally, the combined spectra 
and la error arrays were normalized to the flux density, S v , of 
the quasars. S v was taken as a power-law fit to the flux densities 
presented in Fig.Q 

The combined spectra, normalized by the continuum, are pre- 
sented in Figs. |2jj5] Each contiguous spectral segment is plot- 
ted separately and all regions are plotted with the same linear 
frequency-scale. For clarity, the data have been boxcar-smoothed 
over 3 channels (we analyse the unsmoothed spectra in Section[3}- 
From Fig.|2|it is clear that we only barely detect the continuum of B 
0500+019. However, B 0648-165, B 0727-115 and B 1213-172 
have S/N as high as ~10 per (unsmoothed) 0.7 MHz channel. We 
obtained the largest spectral coverage on B 1213 — 172, a total (dis- 
continuous) bandwidth of 23.2 GHz. 



3 SEARCHING FOR ABSORPTION SYSTEMS 
3.1 Search algorithm 

Table |2| lists the transitions detected in the 4 known mm-band 
molecular absorption systems. From this table we selected a set 
of 'commonly' detected molecules for which to search in our ob- 
served spectra: CO, HCO+, HCN, HNC, CS and CN. We also 
searched for some common isotopomers of these molecules, which 
have lower terrestrial abundances, to provide potentially greater 
redshift space coverage and, as explained below, to help rule out 
some candidate absorption systems: 13 CO, C 18 0, C 17 0, H 13 CO+, 
H 13 CN. For each molecule, we searched for all transitions lying 
b elow 1000 GHz (rest-frame) listed in the molecular line database 
of lPickett et alJil99^) 4 . With this fiducial set of transitions we ob- 
tain complete redshift coverage up to 2 ~ 10 for B 0648—165, 
B 0727—115 and B 1213 — 172. Since the frequency coverage for 
B 0500+019 is significantly smaller, several small redshift gaps 
(i.e. Az < 0.05) appear for 0.19 < z < 3.84. However, the frac- 
tion of redshift space covered over this range is still high, ~85%. 
The redshift coverage is complete for all other redshifts below 
2 « 10. 

In Fig.|6|we illustrate the redshift coverage up to 2 a b s = 3.0 
achieved for B 1213 — 172 using only the strong (i.e. highest terres- 
trial abundance) isotopomers in the fiducial set. From the references 
in Table[2]it is clear that the transitions of HCN, HNC, CS and CN 
typically present a lower optical depth for absorption compared to 
CO and HCO + . Therefore, our estimate of the redshift range cov- 
ered using the fiducial set of transitions may appear optimistic. The 
lower two panels in Fig.[6]compare the redshift coverage when us- 
ing all the strong isotopomers with that obtained using only CO 
and HCO + . Though a noticeable decrease in the redshift coverage 
is clear, just using CO and HCO + as a search diagnostic still allows 
much of the redshift space to be scanned. Note, however, that the 
third panel (from the bottom) in Fig.|6|shows that the redshifts cov- 
ered by CO rarely overlap with those covered by HCO + . Thus, if 
one is to identify absorption systems with such a redshift scanning 
technique, one is biased toward finding systems containing strong 
HCN, HNC and perhaps CS and CN. This bias can only be avoided 
by observing wider-band spectra. 

Using the 1 a error array for each spectrum, we slid a win- 
dow of width iVwin along the spectrum to search for absorption 
features above a significance level of iV CT standard deviations. Note 

4 Available at http://spec.jpl.nasa.gov 



Table 2. Transitions already detected in absorption at high red- 
shift. 



Transition 


Quasar™ 


Reference 6 


CO(0-1) 


B 


WC94b 


CO(l-2) 


A,C 


WC95; WC96a 


CO(2-3) 


A,C 


WC95; WC96a 


CO(3^t) 


D 


WC98 


13 CO(l-2) 


A,C 


WC95; WC96a 


C 18 0(1-2) 


A 


CW95 


HCO+(0-1) 


B,D 


WC97; CM02 


HCO+(l-2) 


A,B,C,D 


WC95; WC96a,b,c 


HCO+(2-3) 


B,C,D 


WC96a; WC96c; WC97 


HCO+(3^1) 


C 


WC96a 


H 13 CO+(0-1) 


D 


CM02 


H 13 CO+(l-2) 


D 


WC98 


HCN(l-2) 


A,B,C,D 


WC95; WC96a,b; WC97 


HCN(2-3) 


B,D 


WC96c; WC97 


HC 3 N(4-5) 


D 


CM02 


H 13 CN(l-2) 


D 


WC98 


HNC(O-l) 


D 


CM02 


HNC(l-2) 


B,C,D 


WC96a,c; WC97 


HNC(2-3) 


B,C,D 


WC96a,c; WC97 


CS(0-1) 


A 


CWN97 


CS(2-3) 


D 


WC96c 


CS(3^1) 


D 


WC96c 


N 2 H+(l-2) 


D 


WC96c 


N 2 H+(2-3) 


D 


WC96c 


H 2 CO(lio-2n) 


D 


WC98 


H 2 0(ortho) 


A 


CW98a 


C 3 H 2 (loi-2i 2 ) 


D 


CM02 


LiH(CM) 


A 


CW98b 



"Quasars: A = TXS 0218+357, B = PKS 1413+135, C 
= TXS 1504+377, D = PK S 1830 -211. ''References: 
CM02 ICarilli & Mentenll2002l). CW95 ICombes & Wiklinl 

S 1995ft. CWN97 ICombes. Wiklind & Nakai Il997t). CW98a 
ICombes & Wiklind Il998al). CW98b ICombes & Wiklind 
Il998blt. WC94b Iwiklind & Combesl 1 1994th WC95 
1 Wikli nd & Combesl Il995l). WC96a Iwiklind & Combel 
Il996ali. WC96b Iwiklind & Combesl Il996th. WC96c 
Iwiklind & Combes! Il996di. WC97 Iwiklind & Combesl 
Il99ii . WC98 Iwiklind & Combeil99£l) . 

that the apparent significance of a candidate feature will be influ- 
enced by any correlations between the flux in adjacent channels. 
We address this problem in the simulations described in Section 
13.21 Once an absorption feature is identified we associate it with a 
particular transition, allowing us to assign to it an absorption red- 
shift, 2 a bs, with an appropriate error, fcabs, which is taken as the 
redshift interval corresponding to the half-width of the feature. The 
'width' is defined here as the number of channels over which the 
feature remains significant at a level ^ JV CT standard deviations. 

To identify an absorption system we search for another absorp- 
tion feature which, when associated with a different transition, has 
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Figure 6. Schematic diagram illustrating the redshift coverage for B 
1213—172 using only the strong (i.e. highest terrestrial abundance) iso- 
topomers in the fiducial set. The upper six panels show the redshifts covered 
by each molecule given the frequencies observed (see Fig.|5j. The verti- 
cal scale shows the lower-level rotational quantum number, j. The third 
panel from the bottom summarizes the coverage from all j'-levels for each 
molecule. The panel marked 'Total' shows the coverage from all molecules 
and the lowest panel includes only the coverage from CO and HCO+ . See 
text for discussion. 

a redshift z a b s ± 5z a bs- For computational convenience we impose 
an arbitrary but conservative redshift upper limit for absorption sys- 
tems, z a bs ^ 5. This upper limit can be extended to higher redshifts 
but our conclusions below are unchanged. Having identified a po- 
tential absorption system we can find the significance of absorption 
for all other transitions in the set defined above which fall within 
the bandpass of the spectrum. 

Some of the candidate absorption systems can be immediately 
rejected using self-consistency conditions imposed by the relative 
'strengths' of all available transitions. For example, if the detec- 
tion of the absorption system relies on significant absorption in 
13 CO(l-2) but there is no significant absorption in CO(l-2) then 
we can rule out the putative detection. To form a rigid selection 
criteria, we found the integrated flux density, Smt, and 1 a error, 
SSint, for the 'host' line [i.e. the weaker but higher terrestrial abun- 
dance line; CO(l-2) in the example above] and the significantly 
'detected' line. In the analysis below we ruled out candidate sys- 
tems if 

StT > St^ + 2 (SShT + 6S&) . (2) 



Table 3. The most significant double-line detections for each 
quasar in our sample. All transitions which fell into the band-pass 
of our observations are given and column 4 shows the signifi- 
cance of each detection (not corrected for correlations between 
data channels; negative values indicate flux values above the con- 
tinuum). The final column gives the width in channels (as defined 
in Section l3~T] . A, of the two significant features responsible for 
the putative detection. 



Quasar 




Transitions 


Sig. (<r) 


A 




2 8 1483 


pen n 1 1 ■» 


l.O 


2 






CS(16-17) 


— 1.0 








C 18 0(2-3) 


1.8 


2 






n\ A ! ^J— O) 


9 1 










0.6 




B 0648—165 


2.02869 


HCO+(4-5) 


4.3 


2 






HCO+(2-3) 


0.8 








1 1 v. yz. — j f 


1 








H 13 CO+(4-5) 


—0.8 








13 CO(3^1) 


0.9 








C 17 0(3-4) 


0.0 








C 18 0(3^1) 


1.9 








HCN(2-3) 


5.4 


5 






HCN(4-5) 


0.0 








HNC(2-3) 


0.3 








CS(4-5) 


0.1 








CS(8-9) 


0.9 




B 0727-115 


1.51374 


CN(l-2) 


4.2 


2 






13 CO(l-2) 


0.1 








C 17 0(l-2) 


3.0 








C 18 0(l-2) 


4.9 


4 






H 13 CO+(3^1) 


0.4 








HCN(3^t) 


-1.3 








HNC(3-4) 


0.8 




B 1213-172 


1.81705 


C 17 0(l-2) 


1.1 


2 






CN(l-2) 


5.0 


4 






CS(4-5) 


2.1 








CS(7-8) 


1.4 





With the target transitions selected above, the fraction of absorption 
systems rejected in this way was typically ~20%. Relaxing the cri- 
teria in equation|2|does not alter our main conclusions below. 

3.2 Search results 

Consistent with a visual inspection of Figs. 2-5, no strong (i.e. 
t ~ 1) absorption systems exist in our data. For all 4 quasars in 
our sample, no single-line absorption features were identified with 
N a ^ 6.5 a for 3 ^ N^in ^ 31. No absorption systems (i.e. 
double-line features) were observed for N a ^ 5.4 a. We list the 
strongest absorption system candidates in Table|5| 

Our search technique can be used to detect weaker absorption 
lines/systems by lowering the rejection limit, N a . However, this 
results in detection of large numbers of single and double-line fea- 
tures. For example, Fig. shows the number of features detected 
in B 1213—172 if we set JV CT = 3.5 a (solid circles). To be sure 
that most of the these 'weak candidates' are spurious (i.e. the re- 
sult of noise) we constructed synthetic spectra with the following 
procedure: 

(i) For each quasar, we modelled each 1 GHz integration as 
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absorption at arbitrary, unknown redshifts. After defining a set of 
commonly detected molecules, we searched for absorption with 
nearly complete redshift coverage up to z a b s = 5. No candidate 
absorption systems (i.e. two independent features attributable to 
transitions with a common redshift) could be identified where both 
features had a significance > 5.4 a. 

Simulations indicate that systems identified with lower signif- 
icance are consistent with noise. However, as a reminder that we 
cannot completely rule out these putative detections, we provide 
the best absorption system candidates in Table|3| These can be re- 
garded as priority targets for follow-up observations. For our high- 
est S/N spectrum, that of B 0727 — 115, the 3 a optical depth limit is 
~0.3 per (unsmoothed) channel. This is comparable to the optical 
depth limits reached in searches for various molecules in damped 
Lyman-a absorption systems (see lCurran et alj20"0 2 and references 
therein) and compares well with the optical depths of many tran- 
sitions in the 4 known high-z mm-band molecular absorbers (see 
references in Table|2j. Therefore, despite our null result, it is clear 
that wide-band correlator systems can be used to efficiently search 
for molecular absorption systems at unknown redshifts. 



Figure 7. Number of single- and double-line matches found for B 
1213—172 (solid circles) compared with the number found in simulated 
data (crosses) for a range of -/V w i n and N a = 3.5 <r. The 1 a error bars are 
derived from the rms of 200 simulations and the crosses represent the mean 
value. 



Gaussian noise with rms per channel equal to that of the real quasar 
data (i.e. the error arrays for the real and synthetic spectra were 
identical). 

(ii) Each synthetic spectrum was convolved with a Gaussian 
with FWHM equal to that of the autocorrelation function of the 
quasar data. 

(iii) We produced a final combined spectrum using the same 
procedure used for the real quasar spectra (see Section l2~2l . 

The effect of step (ii) is to introduce positive correlations between 
the flux density in neighbouring channels. We confirmed that the 
amplitude and range of correlations match those found in the real 
quasar data by comparing the autocorrelation functions for both 
data sets. 

In Fig. Q we compare the number of single- and double-line 
detections in our spectrum ofB 1213—172 with a Monte Carlo 
simulation using synthetic spectra produced with the above proce- 
dure. We note that the number of lines identified in the real and 
synthetic spectra are the same within the standard deviation of 
the simulations. We verified that this is true for the entire range 
3 < iVwin < 31 (i.e. line widths 5kms _1 < Av < 80kms _1 ) 
and for all rejection limits 3.0 a < N a < 5.5 a, suggesting that 
our simple model of the SEST data is adequate. Thus, we conclude 
that we have not detected a large number of weak absorption sys- 
tems in the quasar data, though we cannot rule out single [or small 
numbers (~5) of] weak systems. 



4 DISCUSSION 

We have selected 4 millimetre-loud quasars which have not been 
detected optically (ms > 20), possibly because an intervening 
absorption system causes a large visual extinction. Such dusty sys- 
tems could lead to detectable mm-band absorption and so we have 
performed wide-band millimetre-wave spectral scans to search for 
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